Behavior of quantum entropies in polaronic systems 
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Quantum entropies and state distances are analyzed in polaronic systems with short range (Hol- 
stein model) and long range (Frohlich model) electron-phonon coupling. These quantities are ex- 
tracted by a variational wave function which describes very accurately polaron systems with ar- 
bitrary size in all the relevant parameter regimes. With the use of quantum information tools, 
the crossover region from weak to strong coupling regime can be characterized with high precision. 
Then, the linear entropy is found to be very sensitive to the range of the electron-phonon coupling 
and the adiabatic ratio. Finally, the entanglement entropy is studied as a function of the system size 
pointing out that it not bounded, but scales as the logarithm of the size either for weak electron- 
phonon coupling or for short range interaction. This behavior is ascribed to the peculiar coupling 
induced by the single electron itinerant dynamics on the phonon subsystem. 



CO 



I 

C 

O 

o 



> 

On 

o\ 
o 
o 



INTRODUCTION 

In the last years quantum entanglement has attracted 
enormous interest as a key physical resource at the basis 
of quantum information processing. [l[ In particular, at- 
tention has been devoted to clarify and quantify quantum 
entanglement in many-body systems since entanglement 
measures provide insights on the quantum correlations 
of many-body functions. Q For example, entanglement 
entropy and its scaling properties are currently used in 
order to better characterize quantum phase transitions. 
In the case of local couplings between degrees of freedom, 
the entropy of the reduced state of a subregion grows like 
the boundary area of the subregion, and not like its vol- 
ume, that is known as area law. [3| However, close to 
a transition, this law is violated since the entropy be- 
comes divergent as a function of the system size. The 
size scaling very often follows a logarithm law. 

Tools given by quantum information have been espe- 
cially important for mesoscopic systems where several 
schemes have been proposed for detection and measure- 
ment of entanglement. Q Very recently, a microwave- 
frequency mechanical oscillator has been cooled to its 
ground state with high probability and coupled to a quan- 
tum bit. ||| This coupling preserves the quantum 
states and allows a time-domain control of the system. 
In these experiments the maximum number of phonons 
in the relevant mechanical mode is very low, so that one 
expects that the coupling between mesoscopic resonator 
and quantum bit is not strong. Moreover, a scenario has 
been proposed in order to detect entanglement of a me- 
chanical resonator and a qubit in a nanoelectromechanic 
setup. 

The realization of such devices where many quantum 
bits or multiple electronic states are coupled to meso- 
scopic or macroscopic mechanical objects is far from be- 
ing easy. Therefore, it is important to make theoreti- 
cal studies of such systems elucidating the role of the 



coupling between electrons and oscillators and their en- 
tanglement properties. In particular, the dimension of 
the system becomes a very important parameter for the 
analysis. Instead of considering artificial devices, one can 
analyze compounds already existing in nature where a ar- 
ray of microscopic oscillators is present in mesoscopic or 
macroscopic systems. Therefore, the aim of this work is 
to study from a quantum information perspective a sys- 
tem relevant in many areas of condensed and nanoscopic 
matter: the polaron, i.e. a single electron (with many 
accessible states in a bulk crystal or a quantum dot) in- 
teracting with lattice phonons. 8] It has been proved 
that it does not show any self-trapping phase transition 
with increasing electron-phonon (el — ph) coupling , but 
a crossover between solutions with small extended (large 
polaron) and large localized (small polaron) lattice de- 
formations. [9] Up to now, only quantum entropies have 
been analyzed for polaronic systems with small size and 
short range (SR) el—ph interactions. 13, II | In addition 
to quantum entropies, in this work, distance measures are 
studied in polaronic models with arbitrary size taking the 
free electron as reference state. We will use a variational 
approach that is very accurate in all the coupling and 
adiabaticity regimes. 12, By means of these tools, 



the precise position of the crossover region between weak 
and strong coupling regime is identified. 

Polaron studies extend the analysis of simpler spin- 
boson models. [l4| The entanglement for polaron is very 
interesting also for other reasons: it can be considered as 
a measure of decoherence of the electronic state due to 
the coupling with phonons or dephasing of the phononic 
state induced by the interaction with the single electron. 
From this point of view, studies of polaronic entangle- 
ment are relevant in mesoscopic systems such as quantum 
dots in polar semiconductors that have been proposed as 
systems for quantum processing. [15j | Moreover, in these 
systems, the el—ph interaction is not local. Therefore, in 
this paper, we have for the first time analyzed the effect 
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of long range (LR) el — ph couplings on entanglement 
amount. It is found that quantum entropies are strongly 
dependent on the adiabatic ratio and range of the inter- 
action. The final part of the paper will focus on the size 
scaling of the entanglement entropy. Even if the system 
is not critical, the entropy of polaronic systems is un- 
bounded, and it scales as the logarithm of the size cither 
for weak el — ph coupling or for local interactions. 



MODELS AND VARIATIONAL APPROACH 



state energies, but also effective masses and spectral 
weights calculated with this approach have been com- 
pared with the results of numerical approaches finding 
excellent agreement. The trial wave functions are trans- 
lational invariant Bloch states obtained by taking a su- 
perposition of localized states centered on different lattice 
sites 
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The prototype model with SR local coupling is the well 
known Holstein one .llal while that with LR interactions 
is the Frohlich one. [17| In particular a discrete version of 
the Frohlic model will be studied. [3] The Hamiltonian 
of the polaronic systems is 



H 
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S^(R n ) = J2 [hf{k)^ An + h.c. 



(7) 



In Eq.(fT]) c\ (cj) denotes the electron creation (annihila- 
tion) operator at site i, whose position vector is indicated 
by Ri, and the symbol <> denotes nearest neighbors 
linked through the transfer integral t. The operator a\ 
(a,) represents the creation (annihilation) operator for 
phonon on the site i, and loq is the frequency of the opti- 
cal local phonon modes. Both the SR and LR models can 
be described by the general el — ph Hamiltonian H e i- p h 



H e i-ph = aujQ ^2 f(\Ri - Rj\) c \ci [aj 



(2) 
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where a controls the strength of el — ph coupling, and 
f(\Ri — Rj\) is the interacting force between an electron 
on the site i and an ion displacement on the site j. The 
units are such that the lattice parameter a = 1 and h = 1 . 
The Hamiltonian ([!]) reduces to the Holstein model for 



In the last equations, N is the number of lattice sites (cor- 
responding to the dimensionality of the electron Hilbert 
space), the apex i = w,s indicates the weak and strong 
coupling polaron wave function, respectively, |0 > de- 
notes the electron and phonon vacuum state, tip (k) are 

the phonon distribution functions and 4>^\R m ) are vari- 
ational parameters defining the spatial broadening of the 
electronic wave function. For each function, the varia- 
tional minimization becomes accurate extending the elec- 
tron wave function up to a few neighbors. 

The ground state properties are determined by consid- 
ering as trial state \ipj: > & linear superposition of the 
weak and strong coupling wave functions 
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f(\Ri-Rj\) 



(3) 



while in the LR case [l8| the interaction force is given by 



f(\R i -R j \) = {\R i -R j \ 2 + l 



(4) 



Through the matrix element M^, defined as the lattice 
Fourier transform of auJof(\Ri\), one defines the pola- 
ronic shift E p — YlqMl/wo and the coupling constant 
A = Ep/zt, with z lattice coordination number, that rep- 
resents a natural measure of the strength of the el — ph 
coupling for any range of the interaction. Another im- 
portant parameter of polaronic systems is the adiabatic 
ratio 7 = uio/t. 

We adopt a variational approach previously proposed 
12l Il3| | for the study of systems with variable range 
el — ph interactions and arbitrary size. Not only ground 



where > is the normalized wave function and =< 

4 ^|4^ > is the overlap factor of the two wave func- 
tions. In Eq.® Aj: and B-^ are two additional variational 
parameters which provide the relative weight of the weak 
and strong coupling solutions for any particular value 
of k. In the rest of this work, we will study the one- 
dimensional ground state corresponding to k = k = in 
the physically relevant adiabatic regime 7 < 1. 



RESULTS 

The main quantity extracted from the wave function 
is the phonon-traced electron density operator p e i 



Pel = Tr ph [\ipk=o >< ipk=o\ 



(9) 
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where Tr p h denotes the trace over the phonon degrees of 
freedom. In order to analyze the entanglement between 
electron and phonon, one can use the linear entropy Sl 



S l = l-Tr el [( Pe i) 2 ] 



(10) 



where Tr e i stands for the trace over the electronic de- 
grees of freedom. 19] In Fig. 1, we report the linear 
entropy as a function of the coupling A for SR and LR 
interactions. It is zero for a free electron, then it increases 
with A reaching the saturation value 1 — 1 /N that marks 
the transition to the totally mixed state. For the Hol- 
stein model, Sl reaches the saturation value at a value 
of A slightly larger than unity. In the LR case, Sl in- 
creases due to the larger entanglement between electron 
and phonons. Moreover, the crossover between the weak 
and strong coupling regimes is smoother. 13] All these 



features make the linear entropy a very important quan- 
tity to measure the change of polaron features as function 
of the coupling A. 




FIG. 1. Linear entropy as a function of A for SR and LR 
ranges of interaction at 7 = 0.25 and N = 128. In the inset, 
fidelity and spectral weight as a function of A. The fidelity in 
the fully adiabatic regime is also shown. 

Another important quantity for the analysis of pola- 
ronic systems is the fidelity F 



F = Tr e i [(pfree-el)(Pel)} 



(11) 



where Pfree-ei is the free electron density o per ator ap- 
propriate for periodic boundary conditions. [20| The fi- 
delity considered here is a measure of distance between 
the polaron and the free electron state. In the inset of 
Fig.l, we show the fidelity for LR and SR case. It is 1 for 
A = and decreases with increasing A (for SR it is linear 
in the intermediate regime) up to a value close to zero 
in the strong coupling regime. Therefore, in the max- 
imally entangled state, single electron and phonons are 



so strongly coupled that free electron features have com- 
pletely disappeared. In the LR case, due to the larger 
mixing between degrees of freedom, the fidelity shows a 
marked tendency towards strong coupling features. 

In the inset of Fig.l, the fidelity is compared to the 
ground state spectral weig ht Z = I < Vfc=ol4 = 0|0 > 
| 2 ), with c\ electron creation operator in the momentum 
representation. Z measures how much the quasi-particle 
is different from the free electron (Z = 1). A small value 
of it indicates a strong mixing of electronic and phononic 
degrees of freedom. Z and F (in the same inset) share 
the same behavior as function of A for both SR and LR 
case. 




FIG. 2. Norm distance as function of A for SR and LR el—ph 
coupling at 7 = 0.25 and N = 128. The same quantity in the 
fully adiabatic regime is also shown. In the inset, the trace 
distance as function A for the same values of parameters. 

The spectral weight has been used to distinguish qual- 
itatively a crossover regime (0.1 < Z < 0.9) between 
the quasi-free-electron one (0.9 < Z < 1) and strong 
coupling one (Z < 0.1). [H, [l3[ Our analysis shows 
that a new quantum measure, the norm distance, is able 
to characterize in a quantitative way the crossover re- 
gion. The norm distance is defined as the eigenvalue 
norm of the difference density operator a = p e i— p e i-f ree : 
\\cr\\ n = maxi |j7j|, with r\i eigenvalues of a. 19] In Fig. 
2, the norm distance is reported for SR and LR cou- 
plings. In the first case, it is peaked at a value of A 
slightly smaller than unity. The LR case it is much more 
interesting since the distance shows a maximum at about 
A = 0.5, which corresponds to an intermediate value of 
the fidelity and spectral weight. Finally, it is possible to 
evaluate the trace distance as the trace norm of a/2. [l9[ 
This function shows a behavior very similar to the lin- 
ear entropy. From the comparison of the two distances, 
it emerges that the maximum eigenvector of a is always 
a fraction of the sum of all the others. Therefore, the 
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peak shown in the norm distance is to be ascribed to the 
peculiar structure of the wave-function in this regime. 

It is important to access the difference of the results 
between adiabatic regime and fully adiabatic limit. In 
the limit ujq — ¥ 0, the phonon fields are classical and the 
only relevant el — ph coupling is A. Clearly, quantum 
entropies vanish. However, it is still possible to study 
the fidelity F (shown in the inset of Fig.l). F in this 
regime shows an abrupt change at A = 1 for the SR case. 
This corresponds to the self-trapping transition towards 
a very localized state which breaks the translation in- 
variance. In Fig. 2, we report the norm distance in the 
full adiabatic limit for the SR case. It is strongly larger 
than its corresponding quantity with quantum phonons. 
Actually, in strong coupling, one eigenvalue is close to 
unity, that relative to the electron localized on a single 
site. 

Since quantum entropies are zero in the fully adiabatic 
regime, we analyze the behavior of the linear entropy 
close to this limit. In Fig. 3, we report the linear entropy 
for the SR and LR case as a function of the adiabatic 
ratio 7 for two values of the el — ph constant a. The 
entropies get larger with increasing the adiabatic ratio. 
They are strongly dependent on a and on the range of 
the el — ph interaction. Indeed, for large values of a, the 
parameter A increases very fast as a function of 7. There- 
fore, the linear entropy reaches the saturation value close 
to unity. Actually quantum entropies are very sensitive 
to quantum phonon fluctuations, so that they could be 
used as analyzer of the quantum nature of the oscillators. 
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FIG. 3. Linear entropy as a function of the adiabatic ratio 
7 for several values of the el — ph coupling constant a at 
N = 128. 

In the last part of the paper, we will analyze the en- 
tanglement or von Neumann entropy Svn 

Svn = Tr el [( Pel )ln( Pel )} = S$ N = Sf° n , (12) 



in particular its dependence on the size N . The equiva- 
lence between electron and phonon entropy is due to the 
fact that the starting point is a bipartite pure state. [l9} 
Svn is zero for decoupled electron and phonon degrees of 
freedom, and has the maximum value of ln(N), since ./V 
is the dimension of the smallest subsystem, the electron 
one. 

At hxed size, the linear and von Neumann entropy are 
monotonic functions of each other. Therefore, Svn in- 
creases as a function of A up to a saturation value in the 
strong coupling regime when the state is maximally en- 
tangled. We have found that, in this limit, for both SR 
and LR interactions, Svn reaches the maximum value, so 
that it is not bounded but scales as ln(N). This behavior 
in the strong coupling regime is due, in our opinion, to 
the peculiar "long-range" coupling induced by the single 
electron itinerant dynamics on the phonon subsystem. 




FIG. 4. The quantity £ as function of A at 7 = 0.25. In the 
inset the von Neumann entropy as a function of the system 
size N at fixed value of A. 

This result is even more surprising since it is not lim- 
ited to the strong coupling limit. The scaling propor- 
tional to ln(N) is obtained also in the weak and inter- 
mediate regime. At the second order of the perturbation 
theory in the el — ph coupling, in the limit of large N, 
Svn Nphln(N), where N p h is the perturbative aver- 
age number of excited virtual phonons (linear as function 
of A). We have found by means of an accurate fitting pro- 
cedure that, in all the regimes, Svn scales as £ln(N) + 0, 
where £ (see Fig. 4) and a depend on the el — ph cou- 
pling. In the inset of Fig. 4, we report the von Neumann 
entropies at fixed value of A as function of the size N. 
With increasing the el — ph coupling, £ deviates from the 
linear dependence, and in the crossover regime it curves 
towards the value of unity for strong coupling. Indeed, £ 
behaves as the linear entropy shown in Fig.l. 
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One of the main results of this work is that the entan- 
glement entropy is unbounded and scales as a logarithm 
of the size in all the regimes. This result is valid also for 
LR el — ph interaction. Actually, Sy n cannot increase 
more than ln(N) due to the dimensional constraint of 
the electron Hilbert space. This study allows to quantify 
the amount of the entanglement of realistic systems, for 
example quantum dots in polar semiconductors, where 
LR polaronic effects can be important. In these systems, 
the entanglement can be also related to the measure of 
decoherence of the electronic state induced by the inter- 
actions with phonons. As a result of this study, one can 
estimate that the entanglement entropy due to el — ph 
coupling scales as the logarithm of size with a proportion- 
ality constant that varies in a simple way as a function 
of the parameters. 

The experimental detection of entanglement is difficult 
for the bulk, but feasible for nanostructures. Recently, 
an electronic measurement has been suggested in order 
to detect entanglement between a qubit and an oscillator 
making use of an atomic point contact. Q Furthermore, 
again in the case of the interaction between qubit and 
resonator, time domain control has been used in order 
to controllably create a phonon in the resonator and to 
observe the exchange of this excitation between qubit 
and oscillator. [6] Effects of electric fields have been also 
analyzed for polarons in bulk semiconductors. 2l| The 



pulse induces not only coherent lattice vibrations but also 
velocity drift oscillations of the electron. This is again a 
demonstration of the more complex entanglement behav- 
ior between electron and phonon degrees of freedom that 
has been the focus of this work. 

We acknowledge R. Fazio for a critical reading of part 
of the manuscript. 
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